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Abstract 

We study the properties of thermal radiation emitted by a thin dielectric slab, employing the 
framework of macroscopic fluctuational electrodynamics. Particular emphasis is given to the ana- 
lytical construction of the required dyadic Green's functions. Based on these, general expressions 
are derived for both the system's Poynting vector, describing the intensity of propagating radia- 
tion, and its energy density, containing contributions from non-propagating modes which dominate 
the near-field regime. An extensive discussion is then given for thin metal films. It is shown that 
the radiative intensity is maximized for a certain film thickness, due to Fabry-Perot-like multiple 
reflections inside the film. The dependence of the near-field energy density on the distance from 
the film's surface is governed by an interplay of several length scales, and characterized by different 
exponents in different regimes. In particular, this energy density remains finite even for arbitrarily 
thin films. This unexpected feature is associated with the film's low- frequency surface plasmon 
polariton. Our results also serve as reference for current near-field experiments which search for 
deviations from the macroscopic approach. 

PACS numbers: 44.40. +a, 78.66.-w, 05.40.-a, 41.20. Jb 
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I. INTRODUCTION 



A piece of nonmagnetic, linear and isotropic dielectric material with frequency-dependent 
permittivity e{uj) kept at some finite temperature T generates and emits an electromagnetic 
field, which manifests itself as heat radiation, resulting from thermal and quantum me- 
chanical fluctuations. Besides the thermal far field, near-field phenomena associated with 
nonpropagating modes have recently attracted increasing attention Possibly important 
effects have been revealed, such as the emergence of both ternporal and spatial coherence in 
the near field of planar thermal sources due to surface waves j2 -4|. In addition, the influence 
of these surface waves on radiative heat transfer and dispersion forces at the subwavelength 

n 

scale has been investigated py]. In such studies, one usually considers a simple half-space 
geometry, or two half-spaces separated by a narrow vacuum gap. 

The characteristic length scale for absorption of heat radiation by the material is the skin 
depth 



evaluated at the dominant thermal frequency Uth = 2.821 A;bT//1. We use the notation 
ko = u/c, where c is the velocity of light in vacuum; ej-{uj) = e{uj)/eo is the relative per- 
mittivity, with denoting the permittivity of the vacuum. Thermal radiation generated 
inside the material can reach its surface only if it originates from its outermost layer with 
thickness on the order of ciskin- Hence, if the linear dimensions of a given dielectric are signif- 
icantly larger than the skin depth, the emitted radiation preserves no information about the 
material's geometry. In that case, radiating and non-radiating coniponents of the thermal 
electromagnetic field equal those emitted by a dielectric half-space 5|, |6|. 

In this paper we investigate a seminal example which shows that measurable effects occur 
when the above two length scales become comparable, so that the half-space model becomes 
inadequate. We consider an infinite, planar dielectric slab and study the dependence of 
both its thermal far and near field on its thickness d. When d is large compared to dskin, the 
propagating radiation emitted by the slab is described by the Planck-Kirchhoff radiation 
law [5!]. However, when d is reduced below c?skin, two competing trends arise: On the one 
hand, multiple reflections of radiation inside the slab lead to a Fabry-Perot-like enhancement 
of the field; on the other, the radiating source volume is diminished. We work out the 
consequences of this competition and show that, in particular, the near-field energy density 



2 



close to the surface of a metal film can remain finite even when the thickness of that film 
becomes arbitrarily small, as a result of the emergence of a low-frequency surface plasmon 
polariton. 

We proceed as follows: In Section [II we briefly collect the required elements of fiuc- 
tuational macroscopic electrodynamics 7|, and outline in Sec. IIIII the construction of the 
classical dyadic Green's functions for the dielectric slab; these turn out to be significantly 
more complicated than the more often considered ones for a half-space, or for a vacuum gap 

|. Although properties of propagating thermal radiation may also 
_and numerical codes for investigating thermal radi- 
10(, the detailed analytical discussion of these Green's 
functions is of its own intrinsic value. Green's functions for layered media appear in a variety 
of contexts, such as van der Waals forces in multilayer systems ll|, ll2|, magnetic noise in 
conducting slabs |l3|, light scattering and control of spontaneous emission in planar cavi- 



between two half-spaces 
be obtained by more direct means [9 
ation of layered structures do exist 



ties 
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Il5| . or thermal spin flips in atoms chips [16|, to name but a few. Unfortunately, the 
algebra involved in writing down such Green's functions, though not difficult in principle, is 
vexatingly cumbersome. Here we resort to the formalism outlined in Ref . , based on the 
systematic use of vector wave functions [isl, which combines versatility with transparency, 
and which allows us to treat both far-field and near-field effects on equal footing. With 
the help of the Green's functions we then derive in Sec. IIVI general expressions for both 
the intensity of the slab's thermal radiation field and its energy density, deferring tedious 
mathematical details to the Appendix El Section |V] contains an extensive discussion of these 
results. 

For the sake of definiteness we concentrate on dielectric slabs effectuating a coupling 



between plasma-like electron motion and the photon field 
Drude approach, their permittivity is given by 
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19| . such as metals. Within the 



e{uj) = eo 



(2) 



(1 — iwr, 

where denotes the plasma frequency, and r the relaxation time [20|. Since the dominant 
thermal frequency amounts to Wth ~ 1.1 ■ 10^^ s~^ for T = 300 K, and since typical relaxation 
times for metals are on the order of 10~^^ s, one can satisfy the inequality ur <^ 1 in the 
infrared, thus arriving at the Hagen-Rubens approximation 21 1 



Krf + i^ 



(3) 



for the relative permittivity. For metals with comparatively short relaxation time, such 
as Bismuth (tbi ~ 2.3 ■ 10^^^ s at T = 273 K), this approximation is quite good indeed. 
However, a Drude metallic state can also be achieved with conducting polymers; for instance, 
hexafiuorophosphate doped polypyrrole [PPy(PF6)] yields a plasma frequency in the far 
infrared at about 2 ■ 10^^ s~^, combined with an anomalously long scattering time quantified 
as 3 ■ 10~^^ s in Ref. 22|. Hence, while we use this approximation ([3]) for deriving various 
analytical estimates, we rely on the full Drude permittivity (|2]) in our numerical calculations. 
Besides Bismuth, for which ujpT ^ 4.8 at room temperature, we will also consider more 
typical metals, for which UpT is two orders of magnitude larger. 

We consider in subsection IV Bl the radiative intensity emitted by thin metal films, and 
demonstrate that the opposing trends hinted at above result in an optimum film thickness 
which maximizes that intensity. We then study in subsection IV CI the "evanescent" near- field 
energy density, and establish a somewhat counterintuitive result: While the contribution of 
the TE modes to that density vanishes when the film thickness goes to zero, that of the 
TM modes does not, but remains finite and becomes universal, at least within the scope of 
the simple Drude approach. Finally, we briefly spell out some experimental ramifications in 
Sec.lVD 

II. ELEMENTS OF FLUCTUATIONAL ELECTRODYNAMICS 

As is customary, we consider the macroscopic electric and magnetic fields inside the 



dielectric material, E(r, t] 
some appropriate volume 



and H(r,t), obtained by averaging the microscopic fields over 



23| . so that their small-scale, "atomic" fluctuations are smoothed 
out. Since these fields are described by real numbers, one has 

/ + 00 1 
^E(r,a;)e"'-* 

-^E(r,a;)e-'-* + c.c., (4) 

where c.c. denotes the complex conjugate of the preceding term. A corresponding identity 
holds for H(r, t). Hence, it suffices to restrict the temporal Fourier transforms E(r, w) and 
H(r,a;) to positive frequencies u. 

Q 

Following Rytov and co-workers [7|], we describe the connection between the electromag- 
netic field and its sources by augmenting the dynamical macroscopic Maxwell equations by 



fluctuating current fields. For nonmagnetic materials, characterized by the permeability /io 
of the vacuum, only an "electric" current is required, the frequency components j(r, cu) of 
which are regarded as independent stochastic variables. The resulting equations 

VxE(r, w) = ia;/ioH(r, w) (5) 
VxH(r,a;) = -icje(u;)E(r, w) + j(r, (6) 

then adopt the status of Langevin-type stochastic equations, with j(r, u) playing the role of 
a stochastic force. As in the theory of Brownian motion, the correlation functions of these 
"forces" then are of central importance. As a consequence of the fluctuation-dissipation 



theorem, they acquire the forms 
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25| 



j/3(r',u;') ) = 



{ja{r,u;)j^{r',u')) = 4nu E{u, f3) e"{uj) 6^^6{r - r') 6{uj - u') 



{Ur,uj)jp{r',u')) = 0, (7) 

where angular brackets indicate an ensemble average, overbars denote complex conjugation, 
Greek letters specify vector components, and 

^(-'/5) = ^^ (8) 

is the mean thermal energy of a quantum mechanical harmonic oscillator with frequency u, 
omitting the vacuum contribution. As usual, /3 = (/cb^)"^ is the inverse temperature vari- 
able; e"{uj) is the imaginary part of the material's permittivity e{uj) = e'{u) + ie"{uj). 

Since the stochastic Maxwell equations are linear, there exists a linear relationship be- 
tween the fluctuating sources and the generated fields, which we write as 

E(r, oo) = ioofio J G^(r, r', u) ■ j(r', co) , (9) 
H(r,a;) = ioofio ^dV G^(r, r', a;) ■ j(r', a;) . (10) 

The dyadic kernels G'^(r, r', w) and G^(r,r',a;) are referred to as the classical electric and 
magnetic Green's function, respectively. Once these kernels are known for the given geome- 
try, the correlation functions ([7]) allow one to evaluate bilinear expressions of the fields such 

as 

(EJr,t)Hs(r,t)) = ^ rduu^E(u,p)e%u) fd'r' (g^ -W') + c.c. , (11) 



which enter into the definition of basic observables, such as the Poynting vector of the thermal 
radiation emitted by the material, or its energy density. In the following section, we outline 
the construction of such Green's functions for a dielectric layer of finite thickness. Readers 
not interested in the mathematical details of this construction can proceed to Sec. IIVI 

III. CONSTRUCTION OF DYADIC GREEN'S FUNCTIONS 



The electric dyadic Green's function obeys the inhomogeneous vector wave equation [17 1 
V X V X G^(r, r', u) - eG^{r, r', u) = 16{r - r') , (12) 

where 
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= —e,{u) (13) 

is the square of the wave number inside a material with relative permittivity Ej-lu) = e{uj)/eo. 
Besides, the Green's function also has to satisfy the proper boundary conditions for the 
respective geometry. For nonmagnetic materials, as considered here, the magnetic Green's 
function is then easily obtained from the electric one through the relation 

= X G^ . (14) 

IW/iO 

For constructing Green's functions, one starts from solutions V to the homogeneous vector 
wave equation 

VxVxV-PV = 0. (15) 
Assuming cylindrical symmetry in planes orthogonal to the 2;-axis, the appropriate solutions 



are provided by the vector wave functions 



id 



M,„.(ft) = ( T„:MM|=-|(„^)e, - ^J^S'^h„^)A e- (16) 
\^ p l-cosJ op IsmJ J 

yk dp IsmJ k p LcosJ 

+ ^Jn{Xp){''!"]{nv)eAe''^ . (17) 
sm J J 

Here, Jn{\p) denotes an ordinary Bessel function of order n; the upper (lower) trigonometric 
function goes with the respective upper (lower) sign. The integer n = 0,1,2,3,... is a 
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discrete mode index, whereas the real, wave number-hke index < A < oo is continuous. 
This wave-number index refers to propagation orthogonal to the 2;-axis, and determines the 
wave number h for propagation in ^-direction through the relation 



(18) 



observe that, in contrast to A, this wave number h generally is complex. It appears as the 
argument of the above vector functions, whereas the cylindrical coordinates (p, z) of r 
are suppressed. The functions 'M.±rL\ih) are associated with (a-polarized) TE modes, the 
functions N-i-„a(^) with (7r-polarized) TM modes. Besides solving the homogeneous wave 
equation f|T5l) . these functions f|T6|) and f|T7|) also satisfy the useful identities 



(19) 



It should be pointed out that, besides the vector functions of the M- and N-type, there also 
exists a third type denoted L. These functions obey V x L = and are not needed in the 
present macroscopic approach [igI, Q • 
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FIG. 1: Slab geometry considered in this work. The region — d < z < is occupied by a dielectric 
with permittivity £2- The permittivities for z < —d and z > are given by ei and £3, respectively. 
A source aX z = z' within the slab emits radiation with unit amphtude in both directions. One also 
has to account for waves reflected and transmitted at both interfaces, with ampUtudes as indicated, 
both for TE and TM modes. 
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We now consider the geometry depicted in Fig. [T] The layer between the infinite planes 
z = and z = —d is filled by a dielectric with permittivity 62 = B2{u})] in this layer we have 
stochastic source currents generating the electromagnetic field. For later convenience, we do 
not yet assume at this point that the layer is embedded in a vacuum, but rather that the 
region z < —d be occupied by another dielectric with permittivity Si, while the half-space 



al 



ow us to apply the 



2; > be characterized by still another permittivity 53. This wil 
results obtained here to bulk materials covered by thin coatings {26 . 

If there were no boundary conditions to respect at the interfaces, but the layer —d < z < 
were extended to fill the entire three-dimensional space, the dyadic electric Green's function 
with boundary conditions at 2 = ±00 specifying "outgoing" waves could directly be adapted 
from the literature 17|]: For z > z', one then has 

• POO °° 9 X ( ^ 

Go(r,r',a;) = — / dA — ^ M±„,(/i2) ® M'i„,(-/i2) + N±„,(/.2) ® N;„,(-/i2) , 

2 L J 

(20) 

whereas the signs of the arguments ±/i2 of all four vector functions have to be reversed when 
z < z' . Here and in the following, an unprimed vector function M or N always carries the 
coordinates of the observation point r, whereas a primed function M' or N' refers to the 
source point r'. The symbol ® indicates a dyadic (exterior) product. Finally, M-i-„a(^2) ® 
M'±„;^(-/i2) is a shorthand notation for M+„A(/i2) ® M'+„;,(-/i2) + M_„A(/i2) ® M'_„;,(-/i2). 

For the slab geometry specified in Fig. [H this function Go(r, r', iS) has to be modified such 
that, apart from outgoing boundary conditions at 2 = ±00, also the boundary conditions 
at the two interfaces can be implemented. To this end, we divide the space into four zones: 
Assuming the source to be located within the slab, so that —d < 2;' < 0, we refer to the 
region —00 < z < —d as zone J. Then —d < z < z' ys, zone JJ, while the remaining piece 
2;' < 2; < of the slab becomes zone J//, and the positive half space z > is dubbed 
zone IV . Since the tangential component of the electric field is continuous at the interfaces, 
one requires 

X Gf/,,, = e, X Gf,/,^ ; (21) 

continuity of the magnetic field's tangential components leads to the further boundary con- 
ditions 

X V X Gf/,,, = X V X Gf,/,y . (22) 
These requirements can be met with the following, piecewise ansatz: In zone / one has only 
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radiation transmitted through the left interface, directed towards z = — oo. Specifying still 
unknown transmission amplitudes Tj,^ and Tj^^.j for the TE and TM modes, respectively, 
we therefore write the electric Green's function Gf (r, r', u) with r in this zone in the form 



"^Z^- (23) 

T^j,M±„A(-/il) ® M'±„;,(/l2) + r^MN±„A(-/^l) ® ^'^nX^K 

In zone // there is leftward-directed radiation (with unit amplitude) emitted by the source, 
but also rightward- moving radiation reflected (with amplitudes -R^^; and Rtm) from the left 
interface, together with further, leftward-moving radiation reflected (with amplitudes R^e 
and Rtm) from the right interface: 

'^^> n=0 ^ 

M±„A(-/i2) + R'TE^±nx{h2)^ ® M!^nx{h2) + R^TE^^nxi-h^) ® M'±„^(-/l2) 
+ {^±nx{-h2) + RiTM^±nx{h2)^ ® ^'±nx{h2) + R^M^±nx{-h2) ® N'±„;,(-/l2) | • 

In zone /// one has, mutatis mutandis^ the same dynamics as in zone //, giving 

"^'J n=0 



M±„A(/i2) + i?^sM±„A(-/i2)J ® M'i„;,(-/l2) + i?^^M±„,A(/i2) ® M;„^(/1; 
+ ( N±„A(/i2) + i?TMN±„A(-/i2) ) ® ^'±nx{-h2) + i?rAfN±„A(/i2) ® N;„;,(/12) 



whereas zone IV provides a mirror image of zone J: 

1 / A\^^- ^rifi 



Of course, the above ansatz is characterized in mathematical terms by stating that a suitable 
solution of the homogeneous vector wave equation fllSp has been added to the particular 
solution (pi]) of the inhomogeneous equation (IT^ . 

We are now left with eight unknowns T^'etm ^^'^ ^t^mte^ which match the number of 
boundary conditions provided by Eqs. ( 12T|) and fl22|) . since these apply independently to both 



the TE and the TM modes. We concentrate on the radiation field in zone IV, and therefore 
evaluate the transmission coefficients T^^ and T^m- This procedure is quite cumbersome, 
but elementary, so we immediately proceed to the result: The contribution of the TE modes 
to the radiation emitted by the slab into zone IV is determined by 

2/l2 ' 



T^E^ih) ® M'{-h2) 



-(/ii-/i2)e'''^''M(/i3)®M'(/i2 



while that of the TM modes follows from 



£2 



£1 



(27) 



(28) 



For ease of notation, we henceforth omit the mode indices "±nA" from the vector wave 
functions, and use the determinants 



= {h, + h2){h^ + /i2)e-^'^^'^ - (/ii - h2){h; - h2)e''''^ 



(29) 



and 



Dii 



hi— + /i2 



h3- + h] e^''^^^ 
£3 



hi — 

ei 



£^3 



(30) 



Inserting these results ( !27|) and ( l28l) into the electric Green's function ( !26|) . the radiation field 
existing in the half-space 2; > is completely specified, since the magnetic Green's function 
for that zone is immediately obtained from Eq. (IT^ . keeping in mind the relations (jH]). 



IV. THERMAL RADIATION EMITTED BY A DIELECTRIC LAYER 

As an application of the above formalism we evaluate the intensity of thermal radiation 
emitted by the slab into the half-space 2; > 0, as given by the 2;-component of the Poynting 
vector, 

{Sz{p,^,z)) = e^,3y{E(s{r,t)Hj{r,t)) 

= e,^,^ Hdu u'Eiu, /3)e"{co) [ dV (Gfy ■ G%') + c.c. , (31) 

where e^p-y denotes the Levi-Civita tensor. The calculation is not trivial and involves several 
algebraic manipulations also encountered when computing the field's energy density, so we 
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collect the main steps in Appendix |X1 Because of translational symmetry in planes parallel 
to the interfaces, the result does not depend on the cylindrical coordinates p, ip : 



(S.Xz)) = ^"d. ^ ^"dA Ae--3'= (T, + 71,) 



(32) 



where h'l is the imaginary part of the wave number = h'^ + ih'^, and the dimensionless 
transmission coefficients T± and Ty are given by 

4Re(/i3) 



Til 



\D±\' 



Re(/i2)Ax + Im(/i2)5j 

Re(/l2er2)^|| + Im(/l2^r2)5|| 



(33) 
(34) 



l^lll'kraP 

employing the relative permittivities = ^j{^)/^Oj together with the auxiliary, real quan- 



tities 



All 



+ /i2p(e' 
hi— + /i2 



2h'J.d 



1) + |/ii-/i2|'(l-e 

2 



~2h'J,d\ 



hi 



£2 

Si 



h. 



(35) 



and 



= 2Im(^(/ii + /i2)(/^i-/^2)(e-'"^'2'-l) 



B, = 2Im( ( /.i^ + /.2 U /^i- - ) (e"^^''^'^ 



£2 

VV 



1) 



(36) 



In the case of thick layers, where /igC^ oo, one has 



All 



is "■21 
1 



\Di\^ 







LDii 



— > 



/^s!^ + /^2 



'£3 



2 ' 



5|| 



(37) 



SO that one recovers the well-known coefficients which determine the radiation emitted by a 

4Re(/i3)Re(/i2) 



Til 



\h3 + h2\^ 

4Re(/i3er3)Re(/i2^r2) 



(38) 
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The opposite limiting case of very thin layers, with /igC? — and /igC? — 0, is less obvious. 
After some tedious algebra, one finds 



Til ^ 



\hi + h;\^ 
4Re(/i3er3) 





to 







hi — 



(39) 



with ko = u/c. If both £1(0;) 
vacuum, this gives 



£0 and £3(0;) = so that the layer is surrounded by 



Til -> 



/iorkr2|' + A' 



(40) 



Since according to Eq. ( !T8|) one has /ig = for evanescent modes with A > /cq, only prop- 
agating modes with X < ko contribute to the Poynting vector. It is then a simple matter 
to invoke Eq. f l32p for computing the intensity of thermal radiation emitted by a very thin 
dielectric layer into the vacuum: 



2kr2 



(41) 



As may have been expected, this intensity is proportional to the thickness d of the source 
layer; when that thickness goes to zero, there are no sources left and the intensity vanishes. 
It is also noteworthy that the emitted intensity acquires substantial contributions from 
frequency intervals within which the layer's permittivity is close to zero 19|. 

Besides the intensity, a further quantity of interest is the energy density (u) of the elec- 
tromagnetic field. In contrast to the radiative intensity, this quantity is sensitive also to 
evanescent modes. Focussing on dielectrics facing the vacuum, so that £3(0;) = Eq, we then 
have to evaluate 



2 ' ' 2 

A calculation which largely parallels the one outlined in Appendix |A] eventually leads to 



(42) 



{u{z)) 



2{2TTcy 



dAAe-2^3- 



1 + 



(Tx + T||) 



(43) 



where Sr^ = 1 is understood when computing the transmission coefficients fl33l) and ([3 
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V. THERMAL RADIATION EMITTED BY A THIN METALLIC FILM 



A. Transmission coefficients for propagating and evanescent modes 



In order to clarify the physical significance of the preceding formal results (1521) and 
we now express the key quantities specifying the radiative properties of the dielectric layer, 
the transmission coefficients (!33|) and (!34l) . in terms of the Fresnel amplitude reflection 
coefficients. For radiation coming from a medium with permittivity ei and going into one 
with permittivity 62, these Fresnel coefficients are given by 23 1 

12 ^1 ~ ^2 



for TE modes, and 



,12 



hi + /12 
^ £1 ^ 



(44) 



(45) 



for TM modes. Introduction of these quantities necessitates to distinguish explicitly between 
propagating and evanescent modes. In all of this Section we assume £3(0;) = Eq, and thus 
study thermal radiation emitted into the vacuum. Then propagating modes are characterized 
by A < fco = w/c, evanescent ones by A > /cq. Defining four functions 

/ = (1 - e-2^2^) + |r^Y(e"^''^'°' - e"^''2'^) 
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2e 



Imj r^2( 



1) 



(46) 



where both the Fresnel coefficients r^^ of the left interface and, hence, also the functions / 
and g themselves carry either the label "_L" or "||", then rearranging the r.h.s. of Eqs. ( l33|) 
and ( IMl) with 6^3 = 1, we find their equivalent form 

1 



rppr 



h _ ^12^32 g2ih2d|2 



(1 



„32|2 



)f-2lm{r''')g 



(47) 



for propagating modes, whereas the corresponding expression for evanescent modes with 
7 = -\/A^ — /cg is given by 



2Im(r32)/+ (1 



„32|2 



(48) 



again for both types of polarization. In the case of thick layers, /igci ^ 1, these coefficients 
reduce to 



,32 1 2 



T^^ ^ -Re(/i3)Im(r 

7 



32 ^ 



(49) 
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so that one correctly reobtains the famihar expressions describing thermal radiation emitted 
by a bulk dielectric {2!, [s-l?]. 

However, here we are interested in the radiative properties of dielectric films thinner than 
the skin depth of the material under consideration, so that there are two competing length 
scales. This case requires substantially more care: Provided \h'2d\ ^ 1 and <^ 1, the 
above functions / and g reduce to 

/ ^ 2h'^d {1 + \r'r) 

g -Ah'^dReir^"^) (50) 

in linear approximation. For strongly evanescent modes with A ^ A;o one has /i2 ~ iA, and 
hence /ig ^ 1/^2 1- In the case of propagating modes, a similar hierarchy can be established 
only if one specifies the film's permittivity 62(00). In the following discussion, we restrict 
ourselves to simple metals described by the Drude approach, resulting in the permittivity ([2]), 
or in its approximation ([3]), provided the Hagen-Rubens condition ojt <^ 1 can be met for 
all relevant frequencies. 

Given such a Drude metal, the equality A^ + = /co£r2 requires l/igl <^ /ig for thermal 
frequencies. By means of the approximation f lSU]) . one then finds \g\ f for sufficiently thin 
films, so that contributions proportional to g may be neglected. Hence, we have 

^pr ^ 2^4rf _ 3 

\1 - r^^r^^l - 2h'^d)f ^ ^ 

2h"d 2h' 
|1 -ri2r32(l -2/i^'c/)|' 7^11^^^ ^ ^ 

for thin metallic films. In terms of the dimensionless variable C, = 2h2d, the dependence 
of both propagating and evanescent thermal radiation on the film thickness is therefore 
determined by the function 



^^^^ ]l-a(l-OP 

^ (52) 



[l-a'(l-0F + a"2(l-0 
with a = r^2j-32_ 
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B. Dependence of radiative intensity on film thickness 



From here onwards we focus on metallic films in vacuum, so that r^^ = r^^ = r, and 
first estimate that film thickness which results in maximum far-field heat radiation due 
to propagating modes. We consider only modes directed perpendicular to the interfaces, so 



that A = 0, and 

2 



= ( 'tt^) (53) 



,1 + v^. 

for both types of polarization. The Hagen-Rubens condition ur ^ 1 entails e" ^ \e'^\, or 
l^rl ~ e'^ ^ 1, implying both \r'\ ^ 1 and \r"\ ^ 1. With a = (r' + ir")^ one deduces 
a' ^ ^ 1 and \a"\ <^ 1, so that it suffices to maximize, instead of the function F{C,), its 
approximate version 

^"'(« = e^Jir^ ■ («^) 

Since I a" I <^ 1, the maximum is located at 



5 max 



(55) 



givmg 



4Imy^ 

(56) 



for the optimum film thickness. This characteristic length fl56p coincides exactly with the 
so-called Woltersdorff thickness, which quantifies that thickness of a metallic film which 
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29j. Indeed, 



maximizes its absorption for frequencies u in the Hagen-Rubens regime 
it has been demonstrated experimentally that metal films absorb more infrared radiation 
if they are made thinner; for most metals, maximum absorptance is obtained for films less 



than 10 ^ m thick 



30|. In terms of the skin depth ([T]), which takes the form 



4kin = —\ — (57) 



Up V UJT 



2l|, and using ko = ou/c, the Woltersdorff thickness can be 



for metals in the infrared 
expressed as 

dw = C^skin^O • (58) 
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FIG. 2: Radiative intensity S emitted at temperature T = 300 K by a Drude metal with plasma 
frequency Up = 2.1 ■ 10^^ s~^ and relaxation time r = 2.3 • 10~^^ s, as appropriate for Bismuth, as 
function of the film thickness d (in meters) . Vertical lines indicate the Woltersdorff thickness ()56p 
and the skin depth ([T|), respectively. Data are normalized with respect to the intensity Sbb emitted 
by a black body. Also indicated is the bulk value (horizontal line) and the prediction based on 
Eq. gH). 

The observation that a film of thickness (iw also maximizes its emission reflects the fact that 
absorption balances emission in thermal equilibrium. 

In Fig. |2] we display numerical data for the radiative intensity emitted at temperature 
T = 300 K by a Drude metal with parameters Up = 2.1 ■ 10^^ s~^ and r = 2.3 ■ 10~^^ s cor- 
responding to Bismuth, as function of the film thickness d; in all our numerical calculations 
we employ the model permittivity ([2]) without the Hagen-Rubens approximation. Although 
macroscopic electrodynamics will presumably start to fail for thicknesses below 10~^ m, we 
also plot data for even smaller d, in order to clearly bring out the asymptotic trend. This 
example confirms the picture drawn so far: The intensity emitted by the film almost coin- 
cides with that emitted by the bulk material when d exceeds the skin depth, but increases 
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FIG. 3: Radiative intensity S emitted at temperature T = 300 K by a Drude metal with ujp = 
1.4 • 10^^ and r = 4.0 • 10"^^ s, as corresponding to silver (full line), and another such metal 
with ujp = 2.4 • 10^^ s~^ and r = 0.8 • 10~^^ s, as corresponding to aluminium (dotted). Vertical 
lines mark the corresponding skin depths. 

substantially when d is made smaller, reaching a maximum which exceeds the bulk value by 
a factor of about 3.8 at a thickness predicted neatly by the Woltersdorff formula (!56|) . and 
then starting to decrease. However, the regime of linear decrease described by Eq. f HT]) is 
reached only for unrealistically small d. For metals with more typical Drude parameters at 
room temperature, such as silver {ujp = 1.4 ■ 10^^ s~^ and r = 4.0 ■ 10~^^ s) or aluminium 
{up = 2.4 ■ 10^^ s~^ and r = 0.8 ■ 10"^^ s), the maximum is shifted to even lower d, as wit- 
nessed by Fig. 131 so that only radiative intensity increasing with decreasing film thickness 
might be observable in such cases. It is noteworthy that the intensity generated by a thin 
film can exceed the bulk limit by more than a factor of 10. 
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C. Near-field energy density for thin metal films 



While the dependence of the radiative intensity on the film thickness thus conforms to 
expectation, the dependence of the "evanescent" energy density, at some distance z from 
the film surface, on that thickness is more difficult to oversee. This is related to the fact 
that for strongly evanescent modes with A ^ /cq both Fresnel cofficients ( l44l) and (145|) differ 
strongly: One finds 



(59) 



for A/Zco ^ 1, so that 



/ 9, 4a; 
Re(r^) ^ 1 



but 



giving 



Im(r^) ^ — , (60) 



Re(r 



2^ ~ 

^ 16A4 V CO 



Im(ri) ^ 2ujTRe{rl) (62) 

for X/ko ^ 1 in the Hagen-Rubens regime. Hence, for strongly evanescent TM modes we 
have Re(ry) ^ 1 and Im(ry) ^ 1, which are precisely the propositions which have enabled 
us to reduce the transmission function (152|) to the simpler form (154|) for propagating modes. 
Hence, we can immediately adapt the result obtained in Eq. f lS^ : For evanescent TM modes 
with X/ko ^ 1, the optimum film thickness maximizing the energy density close to the film's 
surface is given by 

which implies 

kodY - (64) 

within the Hagen-Rubens approximation fl6UI) . using /ig ~ ^- In Fig- IHwe show a plot of 
the transmission function (152|) for evanescent TM modes, obtained for Bismuth parameters 
with u kept fixed at the dominant thermal frequency Uth for T = 300 K. Also indicated is 
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the locus of maximizing values in the A-ci-plane, as predicted by the approximation (l63l) : 
obviously this approximation works quite well. 




FIG. 4: Transmission function -F(^) with = 2/i2(i, as defined in Eq. (I52p . for evanescent TM 



modes, so that a 



The permittivity is given by the Drude formula ^ with parameters 



corresponding to Bismuth, as in Fig. [21 the frequency w = 10^^ s~^ is close to the dominant 
thermal frequency for T = 300 K. The locus of the maximizing argument ^max is well described by 
the approximate Eq. ([63|) . 



On the other hand, for evanescent TE modes the approximation fl62|) implies \a"\ <^ \a'\ 
so that the transmission function now is cast into the different form 



with its maximum located at 



This predicts 



[l-a'(l-OF 
a'-l 



Re{r] 



2K 



Re(r2) 



(65) 



(66) 



(67) 
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FIG. 5: Transmission function -F(^) with ^ = 2/12^, as defined in Eq. (j52|) . for evanescent TE modes, 
so that a = r^. The permittivity is given by the Drude formula ^ with parameters corresponding 
to Bismuth, as in Fig. [J] the frequency uj = 10^^ s~^ is close to the dominant thermal frequency 
for T = 300 K. The locus of the maximizing argument ^max is well described by the approximate 
Eq. ()67p . Observe how the scales here differ from those in Fig. HI 

as the optimum thickness for strongly evanescent TE modes, which simplifies to 



in the Hagen-Rubens regime. Again, we display in Fig. |5]a plot of the function (152|) . and 
indicate the forecast of the approximation fl67|) for maximum transmission. Comparison of 
Figs, m and |5] reveals entirely opposite trends followed by both types of modes: Whereas for 
TM modes maximum transmission occurs for comparatively small A, unless kod is excessively 
low, the maximizing A grows with k^d in the case of TE modes. 

These opposite trends leave their imprints in the near-field energy density. The energy 
density {u'^^){z) associated with evanescent modes of either type of polarization at a dis- 
tance z from the film's surface is obtained by integrating the spectral density g°^{u;d,z) 




2 



3 



(68) 



20 



characterizing a film of tfiickness d in vacuum, 

POO 

{u'^){z) = / dug'^{u;d,z) . (69) 
Jo 



According to Eq. we liave 



"'CV 



d, z) = TdA AV^-^ , (70) 

witli T^^ = or T^^ = T|[^ as given by Eq. (gH]), and 7 = - kl. Since Eq. (EH) 
states T*'^ = ^(^(2/13/7) (l + |r|2)lm(r) for thin films, and since rn does not depend on A 
for strongly evanescent modes, the spectral density gYioj^d.z) associated with TM modes 
is given essentially by the integrated product A^e~^'^^F(^). Since the factor A^e~^'^^ has a 
well-developed maximum at Amax = V-^? ^^e spectral density is maximized if the factor -F(^) 
is adapted to that maximum. In view of Eq. fl6^ . it follows that 



d,n^[oj,z) = ^— (71) 

is that film thickness which maximizes the spectral density Q^i'^'i d, z) at a given distance z. 
This reasoning is confirmed in Fig. [6], which shows the full density g\\{uj; d, z), including the 
contribution from propagating modes, for z = 10~^ m and some representative frequencies u, 
again using the example of Bismuth. It is important to observe that this density is dominated 
by low frequencies in the limit of thin films; the maximizing thickness is well captured by the 
approximate Eq. ( |7T]) . In marked contrast, the density g±{uj] d, z) effectuated by TE modes 
is maximized at a thickness which becomes larger with decreasing frequency, as depicted in 

Fig.m 

For calculating the near-field energy density {u'^^){z) we scale all wavenumbers by z, 
obtaining dimensionless quantities such as r/ = Xz. The resulting integral over rj can then 
be evaluated to zeroth order in the small parameter koz, provided the error thus committed 
falls into a range of frequencies where it is suppressed by the Bose-Einstein function E{u, /3). 
This requirement is satisfied if z ^ Ath, where 

is the characteristic thermal wavelength at temperature T. With this proviso, we are led to 

/ = (1 - e'^'^'^^'') + |rp(e"^'''^/^ - e'^'^'^^^) 

9 = 0, (73) 
21 
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FIG. 6: Spectral density (w; d, z) for TM modes with fixed frequency oo originating from a metallic 
film at temperatm'e T = 300 K in vacuum. The density is evaluated at a distance z = 10~^ m from 
the film's surface, as function of its thickness d. The permittivity is given by the Drude formula ([2]) 
with parameters for Bismuth. The frequencies considered are uj = 10^^ s~^, 5 • 10^'^ s^^, 10^^ s~^, 
5 • 10^^ s~^, and 10^^ s~^, in the direction of the arrow. The film thickness resulting in maximum 
density is given approximately by Eq. (ITT]) . 



resulting in the near-field approximation 

2 f^^duj E{uj,P) 



{u-){z) 



Jq u (27r)^ 



, 2 -2v ^Mr) 
^ a?7 ?7 e |]^ _ ^2 g-2»yd/z|2 



1 



-2r]d/z 



-2-qd/z 



— e 



-Arjd/i 



(74) 



When the film thickness is still large compared to the distance from the film, so that z <^ 
this formula simplifies considerably and yields 



Jo ^ (2vr)2 Jo 



(75) 



which coincides with the expression for the energy density close to the surface of an infinitely 
thick layer, i.e., of a bulk material: For distances small compared to the film thickness, the 
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FIG. 7: Spectral density g±{uj; d, z) for TE modes with fixed frequency a; originating from a metallic 
film at temperatm'e T = 300 K in vacuum. The density is evaluated at a distance z = 10~^ m from 
the film's surface, as function of its thickness d. The permittivity is given by the Drude formula ([2]) 
with parameters for Bismuth. The frequencies considered are w = 10^^ s~^, 5 • 10^^ s~^, 10^^ s~^, 
5 • 10^^ s~^, and 10^^ s"^, in the direction of the arrow. 

energy density contains no information about that thickness. For metal films, the above 
result reduces to 

1 1 



{u-){z) 



Az (27rc)2 jQ 
96 c^z h 



[z < d) 



(76) 



for TE modes, and to 



z^ Jo CO (27r)2 + 

{z < d) 



1 1 {kBTf 



24 u^Tz^ 



h 



(77) 



for TM modes. In both these cases, the first expression on the respective r.h.s. is valid in 
general, whereas the second one requires the Hagen-Rubens approximation. Thus, for rather 
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short distances z <^ d the total energy density is dominated by the TM modes, exhibiting 
the famihar 2;~^-divergence known from bulk materials 

In contrast, the other limiting case of distances large compared to the film thickness, 
d ^ -2 ^ Ath, gives rise to a fairly counterintuitive feature. The general expression fl71|) then 
takes the form 

^ ^ z'L u; (2vr)2 J, ' \1 - rHl - 2r^d/ zW 

■ 2r]^{l + |rp) . (78) 

For TE modes, Eq. (16T]) guarantees that the denominator appearing here may be replaced 
by unity. This simplification then gives 

where once again the second approximate equality hinges on the Hagen-Rubens condition. 
Observe that this result differs from the previous Eq. ( 1761) for the reverse situation only by 
the factor d/ z '^1, which appears reasonable: The energy density caused by TE modes at 
a fixed distance z <^ Ath decreases when reducing the film thickness. 

However, the situation is more delicate when dealing with TM modes, for which Eq. f l60|) 
enforces Re(rp !^ 1 and |Im(ry)| ^ 1, so that the integrand in Eq. ( ITHj) acquires a small 
denominator. We still may write 

|1 - rj(l - 2r/d/^)|2 ^ |2r/rf/z - ilm(rj)|2 , (80) 

and now have to make a further distinction: Keeping in mind that the decisive 7] are on the 
order of unity, we may set 

\l-rl{l-2r^d/z)\^^{27^d/z)\ (81) 

provided 2d/ z ^ Im(ry) for all relevant frequencies u < k-QT/h. In view of Eq. (l60l) . this 
requires 

d 2k-BT , , 
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It is interesting to observe that this condition vahdating the approximation ( IHT]) can be cast 
into a particularly suggestive form involving only the Woltersdorff thickness fl56|) and the 
characteristic thermal wavelength (1721) : 

- » ^ ; (83) 

z Ath 

in terms of the maximizing thickness (17T]) . this means nothing but 

d > d^^Jk^Tjfi, z) . (84) 

Given this, we find 

- ^ 2^2 ^^^J^' {dmUkBT/h, Z)'^d'^z) (85) 

for films which are not too thin. This expression differs from its counterpart ( 1771) by the 
factor z/d^ 1, which is noteworthy: In the near-field regime, the energy density caused by 
TM modes increases substantially when reducing the thickness of the film, despite the loss 
of source volume. The fact that this feature emerges only for the TM modes indicates that 
it is related to surface plasmon polaritons. Indeed, when the film thickness d is sufficiently 
small, the plasmons associated with the two surfaces couple, splitting into one resonance 
with a frequency that converges to Up for — )■ 0, and a second resonance with a frequency 



that approaches zero 19|, |32|. This is exemplified in Fig. [HI where we plot the local density 
of states -D||(cu, d, z), defined through the relation {g] 

POO 

{uf{z))= dujE{uj,(3)Dii{uj,d,z) , (86) 
Jo 

for a distance oi z = 10~^ m from Bismuth films of various thicknesses. The increase of the 
energy density {uY){z) with decreasing film thickness found in Eq. (185!) occurs when the 
low-frequency surface plasmon polariton comes into the range of the thermal frequencies. 

Our findings are illustrated in Fig. [9l which shows the individual densities {u±{z)) and 
{u\\{z)) at a distance z = 10~^ m from the surface of a Bismuth film, together with their 
sum, again as functions of the film thickness d. Whereas the density caused by TE modes 
decreases monotonically with decreasing d, the one associated with TM modes actually 
increases, and approaches a finite value for d ^ 0. These different trends obeyed by the two 
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FIG. 8: Local density of states (w, d, 2:) at distance z = W^^ m above a Bismuth film of 
thickness d, specified by the line style. For d = 5 ■ 10~^ m, one observes the usual surface plasmon 
polariton resonance ai uj = Lo^l\f2. For lower d, the resonances associated with the two surfaces 
of the film couple and split. The higher resonance approaches Wp for vanishing film thickness, 
whereas the lower one approaches a; = 0. This latter, low-frequency resonance causes the universal 
behavior expressed in Eq. 



types of modes result in a non-monotonic dependence of the total energy density at z on 
the thickness d. 

The limiting case of very thin films is reached when d/z <^ rfw/Ath, or 



d < d^^^{kBT/h, z) 



(87) 



This limit may be hard to realize in practice, and fall outside the regime of validity of 
macroscopic electrodynamics for most materials, but it is nonetheless of conceptual interest. 
It necessitates to approximate the denominator in the expression (1781) in the form 



l-ri 



d 



27]- 



d 



2r]- 



(88) 
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FIG. 9: Total energy density {u{z)) = {u±{z)) + {u^^ (z)), together with the individual contributions, 
at a distance z = 10"^ m from the surface of a Drude metal film with thickness d. The Drude 
parameters correspond to Bismuth; the temperature is T = 300 K. Data are normalized with 
respect to the energy density ubb of a black body of the same temperature. 



with none of the two terms on the r.h.s. being neghgible against the other for all relevant 
frequencies. Inserting, one is confronted with 

Im(r||) 



,,,,,, 2 f^duj E{io,(3) . 2 _2„ 



d?7?7 



Jq lo (27r)^ 



2 2 r^E(.,p) 



\2rid/ z — i Im(r 
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Arjd/z 
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2 _|_ / 2ujz 
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Rescaling the ?7-integral according to 







2^ 2 



2 , ?/3e-2«^ 
a / dy 



we arrive at 



n'^z'^UpTd Jq 







2ujz 

UpTd 



dy 



(89) 



(90) 



(91) 
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Clearly, here the w-integral is dominated by low frequencies in the small-d- limit, as discussed 
before in the context of Fig. [6l consequently, we are entitled to replace E{u, (3) by k-oT. Then 
interchanging the order of integration, and using 



oo / r> \ 2 



finally results in 



{d < drn,^{kBT/h, z)) . (93) 



kBT 



Remarkably, the energy density associated with evanescent TM modes not only remains finite 
when approaching the (formal) limit of zero thickness, but it also becomes independent of 
the metal's parameters. This universal feature stems from the fact that the low-frequency 
surface plasmon polariton resonance depicted in Fig. [8] converges to zero frequency for all 
Drude materials. Nonetheless, the width of this resonance depends on the relaxation time r, 
so that the scale cimax below which the universal behavior appears is proportional to 1/r. 
As a consequence, metals such as gold and Bismuth, which possess nearly identical plasma 
frequencies but substantially different relaxation times, reveal the above universality for 



rather different film ticknesses. As will be discussed in Ref. [27|], the rise of the evanescent 
energy density expressed by Eqs. ( 1851) and (l93l) also occurs with thin metal films coating a 
polar dielectric, but is lost when a thin metal film covers another metal. 

Figure [TU] shows a doubly logarithmic plot of the dependence of the energy density {u_i{z)) 
associated with TE modes on the distance z from the film, summing up evanescent and 
propagating contributions. Again we take Drude parameters for Bismuth at temperature 
T = 300 K, and consider films of thickness d = 10^^ m and 7 ■ 10~^ m, together with 
an excessively thin model example with d = 10~^° m. In the first two cases, one clearly 
recognizes the crossover from the z~ ^-behavior predicted for 2; <^ d by the approximation (1761) 
to the 2;~^-dependence implied by Eq. (!79l) for larger distances, d z Ath- As indicated 
by the arrow, {u±{z)) decreases at fixed z when d is reduced. 

Figure [TT] displays the corresponding plot for TM modes. Here the data for d = 10"^ m 
and d = 7 ■ 10~^ m exhibit the change of slope from —3, as deduced from Eq. (fTfl) for z <^ d, 
to —2, as required by Eq. ( l85l) for d <^ z Ath- Moreover, for the model case d = 10~^° m 
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FIG. 10: Energy density {u±{z)) associated with TE modes for Drude metal films with Bismuth 
parameters at T = 300 K, as functions of the distance z from the film. Film thicknesses are 
d = m (dashed-dotted), 7 • 10~^ m (dotted), and 10~^*^ m (full line). For fixed z in the near- 
field regime, this density {u±{z)) decreases with decreasing thickness. Straight lines correspond 
to the approximate results (I76p and (j79p for the thickest and thinnest film, respectively. Data are 
normalized with respect to the energy density ubb of a black body of the same temperature. 

one observes another crossover from that slope —2 appearing when d^^y^{k-QT / z) <^ d to 
the universal 2;~'^-behavior found in Eq. f l93p in the opposite limit. 

Finally, we display in Fig. [12] the total energy densities for these examples. As a conse- 
quence of the opposing trends obeyed by the TE and TM modes, there is a well-developed 
range of distances around z ~ 10~^ m where the total density varies non-monotonically with 
the film thickness. 
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FIG. 11: Energy density {u^^{z)) associated with TM modes for Drude metal films with Bismuth 
parameters at T = 300 K, as functions of the distance z from the film. Line symbols are as in Fig.llOi 
For fixed z, the density (n||(z)) increases with decreasing thickness. Straight lines correspond to 
the approximate results ()77p and ()85p for the thickest and thinnest film, respectively. 

VI. CONCLUSIONS 

A theoretical treatment of heat radiation and thermal near fields generated by thin dielec- 
tric slabs within the framework of macroscopic fluctuational electrodynamics hinges on two 
basic ingredients. On the one hand, Maxwell's equations have to be solved with the bound- 
ary conditions imposed by the slab geometry; on the other, the dielectric permittivity e{uj) 
for the slab's material is required. We have addressed the first problem by constructing the 
dyadic Green's functions for the slab in Sec. Illlt and stated general expressions for both the 
intensity of the heat radiation and its energy density, valid for any prescribed permittivity 
e{uj), in Sec.IVl 

The restricion to metallic films, with dielectric response solely due to free particle-like 
electron motion, leads to an intricate competition of length scales. Besides the geometrical 
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FIG. 12: Total energy density {u{z)) = {u±{z)) + (m||(z)), for the same example cases as studied 
in Figs. [10] and [TTJ Observe that the energy density varies non-monotonically with film thickness 
for distances around z ~ 10~^ m. 

thickness d of such a film, and the skin depth ([T]), the Woltersdorff thickness (iw defined in 
Eq. fl56|) and the characteristic thermal wavelength Ath come into play. As a consequence, 
the dependence of the thermal near-field energy density on the distance z from the film's sur- 
face is not characterized by a single exponent, but different exponents dominate in different 
regimes. In particular, we have shown that for d/z <^ ^wZ-^th the energy density is domi- 
nated by the universal contribution ( 193|) brought about by TM modes. Our conclusions are 
based on the Drude model (E]) for the permittivity, regarding the plasma frequency Up and 
the relaxation time r as fixed. This does not hold exactly for electrical transport in ultrathin 
metallic films, since the resistivity of such a film increases with decreasing thickness 33| . 
but we expect the main features of our discussion to persist at least qualitatively. 

With a view towards laboratory experiments, one might be interested not in metallic films 
in vacuum, but rather in thin metal coatings on a planar surface of another material. This 
case has been anticipated in our general formulae (1471) and ( 148|) . where the Fresnel reflection 
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coefficients for the left interface differ 
account for the bulk contribution 26 



Tom those for the right one, but one still has to 



27|. Qualitatively, however, the Fabry-Perot-like 



effect brought about by thin metal coatings bears interesting possibilities of manipulating 
the intensity of thermal radiation, that is, the amount of energy transported per time from 
a hot body into the vacuum. 

The insights obtained in this work are especially pertinent with regard to recent develop- 



memts in scanning thermal microscopy 



35| . It has been argued that the tip of a thermal 



microscope operating in ultrahigh vacuum is sensitive to the near-field energy density above 



the sample; experiments along this direction are presently underway [36|. Near-field sig- 
natures offered by thin films, or by bulk materials coated with films of varying, exactly 
specified thickness even down to that of atomic monolayers, and their actual observation in 
a suitable experiment could be of value for developing scanning thermal microscopy into a 
quantitative tool for materials science. 

While we have formally extrapolated the realm of macroscopic electrodynamics down to 
even unrealistically small scales, it goes without saying that at some scale deviations from 
the macroscopic description will show up; such deviations now have come within the scope of 



modern experimental set-ups 



361] . In addition, one has to expect important corrections due 



to effects of non-local optical response 
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371 ] . The pursuit of the question at precisely what 



length scale, and in precisely what manner, these deviations start to manifest themselves is 
a quite important task in nanoscale thermal engineering. In this quest, the present study 
may serve as a reference. 
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Appendix A: Calculation of the Poynting vector 

In this Appendix we provide technical details required for the derivation of the expres- 
sion fl32|) for the intensity of heat radiation emitted by a dielectric layer. Computing the 
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magnetic Green's function Gfy according to Eq. (HM from the electric Green's function (126|) . 
and defining the convenient symbols 



hi + hj 

e 



+ hi 



b^l := hi — hj 
W := hi^-^ - h 



1 1 



one starts from the product of the two dyadics and integrates, obtaining 

(2 - 5„,o)(2 - 5n'fl) 
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oo roo 



n,n'=0 



XX' 



+ 



■ M(/i3) ® N(/i, 



3) ■ 



la 



N(/l3)®M(/l3)- 



N 



'N 



l2\2^-2hi'd T++ 
e J]vf 



Here we have introduced the source integrals 
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where the integration extends over the dielectric layer labeled "2", that is, over the slab 
—d < z <Q occupied by the sources (cf. Fig. [1]); the symbol A^'^, is given by 



(A4) 
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Likewise, one has 
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The emergence of the 8 different source integrals flA3p and flASp is a characteristic comph- 
cation brought about by the finite thickness of the slab; it can be seen as a consequence 
of multiple reflections inside the slab. Inserting the right hand sides of these integrals, the 
longish expression (lA2p is reduced to 



2-6. 
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In the next step, we evaluate the vector product appearing in the definition of the Poynting 
vector ( 13T|) . Because of translational invariance in the planes orthogonal to the 2-axis, it 
suffices to consider p = only. With the help of 



(M(/l3) ® N(/l3)) 



p=0 



^3A^ 

2k. 



and 



2k. 



(N(/Z3)®M(/13)). 

and exploiting the elementary but important relations 
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where ko = u/c, one is led to 
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Pi 
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47rcu/io Jo 
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(AlO) 
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with the symbols A±, A\\, B±, and as specified in Eqs. (l35l) and (136|) . Inserting this 
expression (lAlOp into Eq. f l3T|) . adding the complex conjugate, and utilizing the identity 

= ^ (All) 

for the imaginary part of the permittivity, one finally ends up with the result stated in 
Eq. (IHID. 
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